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We derive a well-behaved nonlinear extension of the non-relativistic Liouville-von Neumann dy- 
namics driven by maximal entropy production with conservation of energy and probability. The 
pure state limit reduces to the usual Schroedinger evolution, while mixtures evolve towards maxi- 
mum entropy equilibrium states with canonical-like probability distributions on energy eigenstates. 
The linear, near-equilibrium limit is found to amount to an essentially exponential relaxation to 
thermal equilibrium; a few elementary examples are given. In addition, the modified dynamics is 
invariant under the time-independent symmetry group of the hamiltonian, and also invariant under 
the special Galilei group provided the conservation of total momentum is accounted for as well. 
Similar extensions can be generated for, e.g., nonextensive systems better described by a Tsallis 
q-entropy. 



I. INTRODUCTION 



^ , A number of recent, independent experiments have provided impressive bounds on the possible deviations from 
a linear and unitary propagation of pure quantum states, at least on a laboratory accessible space-time scale. The 
limits imposed in this way on potential generahzations of the standard unitary quantum equations of motions, as 
sought in relation to Hawking's blackhole evaporation process , are likewise severe. Certainly, there always remains 
the possibility of modified dynamical laws on the (inaccessible) Plank scale |^ , as well as under the extreme physical 
. environment characteristic of singular cosmological phenomena. Related models of open system dynamics due to 
' alleged statistical perturbations, e.g. from the space-time foam, have enjoyed considerable attention lately But 
, in case the unitarity of pure state propagation holds under universal conditions, one is necessarily lead to a quest for 
O • genuine nonlinear extensions for isolated systems, possibly involving an explicit arrow of time. Indeed, it has been 
^ I ^ pointed out in a fairly general ansatz |^,^ that if the pure states happen to be attractors of a nonlinear evolution, 
^ then testing the unitary propagation of pure states alone cannot rule out a nonlinear propagation of mixtures. This 
situation has been noted recently in the context of certain nonlinear Lie-Poisson dynamics 1 1] , wherein pure states still 
^ , propagate in the usual hamiltonian way, while density matrices evolve nonlinearly, but preserving a time-independent 
spectrum. Unfortunately, the underlying physics remains rather obscure in these theories and the selection of particular 
realizations relevant to various experimental setups is, in general, a matter of guesswork. 

In the following we show that a physically meaningful nonlinear extension emerges when the fundamental postulates 
^ . of quantum mechanics are supplemented by the first and second principles of thermodynamics, at the sole expense 
' of ignoring the constraint of a linear, unitary evolution in time. The result is a largely irreversible, highly nonlinear 
generalization of the non-relativistic quantum Liouville equation, of a form closely related to the ansatz of Ref. [^(but 
not in the Lie-Poisson class), which features a number of rather intriguing properties. In particular, pure states are 
still propagated unitarily into pure states according to the usual (time-reversible) hamiltonian dynamics. The same is 
true of mixed states characterized by an initial equiprobable distribution on a (finite) set of uncorrelated (orthogonal) 
states. Non-pure states evolve so as to maximize the entropy production at each moment in time and to reach station- 
ary states of maximum entropy (or minimum entropy production, according to Prigogine's nonequilibrium principle 
0) on the shortest path in the appropriate state space. Precisely, mixed states arbitrarily distributed on a finite set 
of uncorrelated states evolve into mixed states distributed on an equal number of uncorrelated states, have a time- 
dependent eigenspectrum and eventually attain stationarity on a subset of energy eigenstates. A similar statement can 
be inferred, by extension, for mixtures of an infinite set of uncorrelated pure states. It follows as well that the proba- 
bility distribution at equilibrium, on (a subset of) energy eigenstates, has a canonical-like dependence on the energy 
eigenvalues. For mixtures with an infinite energy range, the corresponding temperature is, of course, strictly positive. 
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whereas for mixtures of a finite set of pure states the stationary state may display a "negative-temperature" distribu- 
tion, in analogy to systems with a finite-dimensional state space. The above mentioned properties are endorsed by the 
positivity of the underlying evolution equation, which ensues by construction despite the high degree of nonlinearity 
involved. The nature of this essentially irreversible propagation becomes evident in the close-to-equilibrium limit, 
when the matrix elements of the density operator between energy eigenstates are found to undergo simple exponential 
decays to the canonical equilibrium values. Finally, proper (non-relativistic) invariance and conservation properties 
under the symmetry group of the hamiltonian are also accounted for. However, in the absence of an explicit general 
law of entropy increase, the time scale for thermal relaxation is set by one multiplication factor, a scalar functional, 
which is yet to be given a specific expression. 

Unlike the nonlinear Lie-Poisson dynamics j^, our framework apparently challenges the notion of separability of 
isolated, non-interacting systems, lack of which has long been thought to be unacceptable Q. We argue, nevertheless, 
that in a nonlinear theory it is necessary to refine the operational definition of isolation and to acknowledge that the 
mutual isolation of two non-interacting systems prohibits entanglement, if individual time-translation invariance is to 
be preserved. When this restriction is properly taken into account in the formulation of the corresponding equation 
of motion, separability can be easily recovered. On the other hand, the case where non-interacting subsystems are 
allowed to develop correlations spontaneously and eventually exchange energy (heat) is shown to correspond in our 
ansatz to the phenomenon of ideal thermal contact. From a precise technical perspective, the effect has its origin in 
that the second principle applies, as usual, to the total entropy of a compound system and not to the entropies of 
individual subsystems. This necessarily results in such redistribution of probabilities and energy as to maximize the 
overall entropy. In physical terms, an ideal gas is allowed to relax spontaneously to thermal equilibrium. 

The formalism can be adapted straightforwardly to cover nonstandard forms for the entropy and energy functionals. 
As immediate examples, we construct a generalization of the Lie-Poisson dynamics with maximal entropy production 
and a nonlinear extension of the standard von Neumann evolution with maximal increase of the nonextensive Tsallis 
q-entropy ^. 



II. THE MODIFIED EQUATION OF MOTION 

Following an earlier suggestion jl^] , the state of a quantum system will be represented by a generalized "square-root" 
7 of the density matrix p, defined by 

P = 77^ • (1) 

In analogy to the common terminology, the operator 7 (not necessarily hermitian) will be called here a state operator. 
Note that the above decomposition is always well-defined, although not unique, for any hermitian and positive definite 
p. On the other hand, to any given 7 there corresponds a unique hermitian and positively defined p. We also adopt 
the standard inner product on the associated Hilbert space of operators, 

(/3|7) - rr(/3+7) , (2) 

such that for 7 normalized, (7I7) — Tr('-f~^^) = I, the average of an observable O becomes the bilinear form 

(7IOI7) = (7+O7) = TriOp) , (3) 

with O the super-operator defined by O, 

0|7) = IO7) . (4) 

It is further convenient to define the tilde-conjugate A of an arbitrary, and not necessarily linear, super-operator 
A 0,by 

(A|a))+ = A|a+) . (5) 

It can be immediately verified that the super-operator A maps hermitian operators a = a+ into hermitian operators 
/3 = /3+ = A\a) if and only if it is tilde-symmetric, A = A. For a super-operator generated by a linear operator, such 
as in Eq (Q) above, the tilde-conjugate is given by 

A|a) = \aA+) . (6) 

In particular, for the hermitian observable O it reads 
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6\a) = \aO) . (7) 

The tilde operation is distributive against the addition and multiplication of super-operators, (A+B) = A + B, 
AB = AB, and is anti-linear against multiplication by scalars, (aA) = a* A.. 

Let us consider now a massive isolated system characterized by an energy operator (hamiltonian) H and a state 

operator 7 (density matrix p = 77^), in an incrtial reference system where its ccnter-of-mass is at rest. We wish to 
find an equation of motion for this system which is first-order differential in time and such that : 
1) Probability is conserved: 

^(7|7) = |tKp)=0 (8) 



or 



(7I7) + (7I7) = 0, 7 = ^7 • (8a) 

2) Energy is conserved (first principle of thermodynamics): 

|(7|H|7) = |TKi/p) = (9) 

or 

(7IHI7) + (7IHI7) = . (9a) 

3) The entropy production is always positive (second principle of thermodynamics in non-equilibrium form), 

|s(t)>0. (10) 

or 

S(t) = -[Tr{p\np) + Trip)] = -[( 7 I ln(77+) I 7) + ( 7l ln(77+) I 7) + 

+ (7l7) + (7l7)] >0. (10a) 

where we adopt the standard entropy expression for a normalized state {Tr(p) = 1) 

S(t) = -ksTr [p(t) Inp(t)] = -ks( 7(t) | In [7(t)7+(t)] | 7(t) ) , (11) 

with ks the Boltzmann constant. 

In order to construct the desired equation of motion, we find it convenient to consider a stronger form of the second 
principle, by requiring that the entropy, as a functional of 7, increase in time along a path of maximum ascent. In 
other words, let the entropy production (10a) be maximized, for any given state 7, against variations of the time 
derivative 7, under the constraints (8a) of conservation of probability and (9a) of conservation of energy. Note that 
the variation of 7 must avoid the simple multiplication by a positive scalar, i.e. a trivial norm increase, since S(t) 
increases then unconditionally. Hence the entropy production must be maximized against the "direction" of 7, that 
is, against derivatives 7 of equal, but otherwise arbitrary norm. This amounts to deriving the equation of motion 
from the following variational principle with constraints 

5{{^\ ln(77+) I7) + (7 I ln(77+) I7) + 2C(7|H|7) + 2r(7|H|7) + 

+ U(7l7) + (7l7)]+-(7l7)} = 0. (12) 

a 

The variation refers to 7 and 7+ only and the form of the Lagrange multipliers ^, a has been chosen for later 
convenience, a and ^ are real scalars on account of their corresponding real functionals, while ( is allowed to span 
complex values. Upon taking the variation of 7, 7+, one is left with 



7) 



^[ln(77+)]|7) + CH|7) + ||7) 



(13) 
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and the hermitian conjugate. Using Eq.(|13D into conditions (gaj) and ( |9aD immediately gives 

1 (7|H ln(77+) |7) + E^ 
ReC = -- = (14a) 



kB(7l7) 



2 i?eC E , 



(14b) 



with (7I7) = 1 , E = (7|H|7)/(7|7) the average energy of th e sy stem, S > the entropy and Ai?^ — (-/\H^\j) — 
the squared energy deviation. One can also check condition (10a) and find that 



(15) 



|0)=ln(77+)|7) + 2CH|7)+ ^|7), (16) 



hence inequality ( |lOa|) is satisfied provided 

cr > . (17) 

In deriving expression ( |l5| ) we used the fact that for \9) as in Eq.(^, and i?e^, £, given by eqs.(p^, it is also true 
that 

(7|H|0) = O, (7|e)-0. (18) 



Let us stress at once that, unlike the usual stationary action principle, our variational principle Eq.(|l2|) does not 
involve variations of functionals over an extended interval of time, but only variations against 7 which are local in 
time, at each given instant t. As a result, the Lagrange parameters C, need only be constants against these same 
variations of 7 and not constants of time or 7 itself. Likewise, condition ( p^ ) for a only guarantees the positivity of 
the entropy production, but does not make S independent of time. Hence all parameters in the equation of motion 
( p^ ) for 7, as well as the entropy production and the entropy itself, are time dependent through their dependence on 7. 
Furthermore, note that Re^, ^ are really functionals of p and H only and therefore are invariant under transformations 
of the kind 

7 ^ 7C/ , UU+ = U+U = I , (19) 

which leave the density matrix unchanged, 

p — > U+-i+ = 77+ . (20) 



Eq.(|13D will be invariant in its entirety under transformation s (|20| ) provided a and ImC^ are likewise invariant as 

functionals of p and H. In that case the entropy production Eq. (|15|) will also be invariant under transformations (|2^) , 
as should be expected on physical grounds. 

Now let us introduce the equivalent equation of motion for the density matrix, starting from 

p = 77+ + 77+ . (21) 

It follows at once that 

p = -a [ p \np + ReC {H - p} - p Tr{p\np)] + i a{ImC)[ p, H] , (22) 



where {,} denotes the anticommutator, as usual. The commutator on the right hand side of Eq.(|22|) obviously provides 
the unitary hamiltonian limit, and the standard Liouville equation suggests 

a(/mC) = i . (23) 

Setting now, for simplicity, Re(^ — > C, the final form of our equation of motion for the density matrix is found to be, 
in common notation, 
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p\np + aP,H-^){H-¥., p}-p 



Tr(p In p) 



where 



Trip) 



(24) 



Cip.H-E) 



1 Tr[{H -E)p\iip] 
'2 Tr[{H - E)2p] ' 



a{p,H ~ E) > 



Tr{p) = const. {= 1) , 



E = 



Tr{Hp) 
Trip) 



= const. 



S = -ks -Tr(p In p) > . 

The scale setting parameter a remains unspecified so far, and will be regarded in the following as a functional of p and 
H. In order to secure that Eq.(p^) is invariant under a scaling p —>■ ap, it must be assumed that cr(ap, H) = a{p, H), 
in which case scaling invariance is verified straightforwardly. Moreover, since eq.(|^) should not show a dependence 
on the zero-point of the energy, it may be assumed also, as above, that a = a{p,H — E). For simplicity, it will be 
understood throughout the following that Tr{p) — 1. 

It is interesting to note that Eq.(p^ can be recovered from a modified form of the nonlinear ansatz proposed in 
Ref.[5], 



a 
T 



fip) - P 



Trjfjp)) 
Trip) 



with the obvious substitutions 



^ ^ cr, f{p) 



p\np + C{H, p} . 



III. FUNDAMENTAL PROPERTIES OF THE NONLINEAR EVOLUTION 

Eq.(p^ secures the hcrmiticity and positivity of the density matrix by construction, since it has been generated 
from an equation for the state operator 7. Conversely, Eq.(p4[) can be easily decomposed into the corresponding 
equations for 7 and 7+ by using the substitution p — 77^, hence the equations of motion for p and 7 are indeed 
equivalent. 

Assuming again a well-behaved a, Eq.(|24|) is seen to be covariant under time- independent unitary transformations. 



p = U+pU 



H 



H = U+HU 



and, in particular, invariant under the (time-independent) symmetry group of the hamiltonian, [U, H] — 0. But an 
observable which commutes with H, [H, O] — 0, is not, in general, an integral of motion. More details on the 
problem follow in Sec. 4. 

It is convenient to absorb the hamiltonian commutator term by setting, in analogy to the usual Heisenberg repre- 
sentation, 



p(t) = exp 



p{t) exp 



Upon substituting expression (|2^), eq.(gj) becomes 

P = -cr [plnp + C {H - E, p} ~ pTr{p\np)] 



(25) 



(26) 
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Now note that for p corresponding to a pure state, p = = \^){'^\, the entropy operator vanishes together with 
the coefficient C, i-e. plnp 0, C,{H — E) ^ 0, such that p(t) = and p(t) = p(0) = |\E')(\1/|, if a is also finite in this 
hmit. From Eq.jpSj) it follows then that a pure state evolves into a pure state according to the usual hamiltonian law: 



p{i) = p^(t) = exp 



I*) (5' I exp 



(27) 



Another situation where the nonlinear evolution reduces to the hamiltonian law is found for uniform (equiprobable) 
distributions Punif, when the eigenvalues of the density matrix are all identical. In this case one has the identity 
Punif ^TLpunif = PunifTr{punif lnp„„i/) and C{Pumf,H - E) ^ 0, wherefrom Punifi^) = 0, Punif{t) = Punif{0) and 



P«m/(t) = exp 



-Ht 



Punif 

(0) exp 



h 



Hi 



(28) 



Recall that under unitary propagation the cardinality of the set of nonzero eigenvalues of the density matrix is 
preserved in time. The same holds true if the density matrix evolves according to Eq.(|2^). In order to see this, 
let Pv = \<l)y){(f)u\ be the projector on some eigenstate of p(t), d - Pi, — pi^Pu, where p^ — Tr{Pi,p) denotes the 
corresponding eigenvalue. Since Tr{p • Pu) = Pv^ multiplying Eq.(|2g) by P^ and taking the trace yields 



Pv^ -<y\Pv^'^Pv + Q:j.(p, i/) p^] 



(29a) 



a.(p, B) = 2 C(p, B) Tr\ P,.{t){H _ E) ] + ^ . 

KB 



(29b) 



Taking pi, In — > for p,y = gives pi/ — and piy(t) — 0, i.e. a zero eigenvalue evolves into a zero eigenvalue. 

As an immediate corollary, density matrices with a finite number of "occupied" state vectors (i.e. a finite number of 
nonzero eigenvalues) are necessarily driven towards a stationary state with a thermal-like distribution on a finite set 
of energy eigenstates. Indeed, in this case the entropy, as a functional of the eigenvalues and under the constraint 
of conserved energy and probability, has a finite absolute maximum. For this reason, and because S(t) > at all 
times, it can only evolve towards a stationary value less or equal to that maximum. But, as will be shown, S{t) — 
implies in fact p = and [p, H] — 0, and the stationary version of Eq.(p3a) gives then the thermal-like distribution. 



Let us prove now that ^(t) — implies stationarity. We begin by making a change of variables, p^ — e > 0, 

such as to write 



_S_ 

k^ 



and 



_S_ 

k^ 



since ^ 1]^ e"^'' = ^ Pu = Also, Eqs.(^) give 



which taken into Eq.(^ produces 



(30) 



(31) 



(32) 



(33) 



Further, use of the explicit expression for ai,, Eq.(p9b), will show that 



2 C VuTr{P„{H - E)) c-"- + —77, e""'^ 

ke 



E 



4C'[rr(P,(H-E))]2 



AC^Tr{P,AH 

KB 



E)) e-'" + - 
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= -2CTr((iJ-E)plnp)+ ( ^ ) - 4 ('Tr {{H ~ Ef py 



S _ / S ' ^ 



-AC:^Tr{{H-E)p)-i^—] =0, (34) 



where we have used the exphcit expression of C, Eq.(|l4|a). Accounting for Eq.(|3j) in Eg. (|33[) shows that 

wherefrom it follows that S = if and only \i r]^ =0 or, equivalently, pi, = 0. Consider now that the system is 
evolving in an asymptotic region where S(t) — s- for all t > . Since necessarily ^ 0, p must be driven by a 
unitary evolution, p(t > to) = U(t)p(to)U+(t). But for 0, Eq.(|l]) gives \T^Pl, = — a,y, which in turn shows that 

plnp = - V p,a,P, = _2 C V p,P^Tr{P,{}l - E)) - V p,P,^ = 

= -2 C{Hd ~E,p}-^p, 

(36) 

where Hd = '^PyTr{PyH) is the diagonal part of H in the eigenbasis of p, [i?D, p] — 0. Introducing the above 
result into Eq.(p6[), one is lead to 

p= -o- C{Hnd, p} , (37) 



with Hnu = H — Hd the non-diagonal part of H relative to p. But Eq.(p7D cannot generate a unitary evolution 
unless Hnd — 0, which implies that stationary entropy over an extended period of time is equivalent to 

Hd^H , (38) 

hence [p,H] = and p — Q. In other words, the density matrix of the system (see also Eq.(p5|)) is stationary and 
also diagonal over energy eigenstates. The explicit form of the occupation probability corresponding to an (occupied) 



energy state of energy E^ follows from Eqs.(29) 



p7 = exp 



-2C«(E,-E)- — 



(39) 



and can be brought to the recognizable thermal form 



p7 = ^e-^E-- , (40) 

with /? = 2^'^'' and Z = — /3E + (S'^'^/ks). Surprisingly, the parameter ( is seen to become at equilibrium, up to a 
factor of 2, the reciprocal temperature f3 — l/ksT. It should be noted, nevertheless, that according to our initial 
assumptions Eq. (^) applies only to a finite number of energy eigenstates and therefore does not refer to a canonical 
equilibrium distribution. More precisely, the sign of i^*^^, and of the generalized temperature T, is not necessarily 
positive. For instance, let the occupied energy eigenstates be labeled by v in order of their increasing energy E^, and 
let their total number be N. If the conserved average energy E is such that 



N 

I 

E > 



1 ^ 



N 



a simple calculation will verify that the entropy will have an (absolute) maximum, corresponding to the equilibrium 
state, on a distribution characterized by a negative hence a "negative temperature". 

At this point, let us examine more closely the restrictive assumption of a finite number of non- vanishing eigenvalues 
for the density matrix. It can be noted that it has entered the argument developed above solely by way of the related 
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assumption of a finite absolute maximum for the entropy, at the given value E for the average energy. However, there 
is good reason to assume that such an absolute maximum exists at least for a large class of distributions over infinite 
sets of (orthogonal) state vectors. If we can extend this "finite absolute maximum" conjecture to all distributions 
with a finite average energy, it becomes possible to generalize the results in Eqs.( ^8|^0| ) and state that the nonlinear 
dynamics described by Eqs.( p^ , p6| ) drives the system towards an equilibrium state on energy eigenstates, with thermal- 
like occupation probabilities. Of course, when the range of occupied energy eigenvalues extends to infinity, relation 
( p] ) can no longer be satisfied for any finite E, and the corresponding temperature can only be positive. 

Finally, we wish to clarify the consistency of the present nonlinear dynamics, which follows a path of maximal entropy 
production, with Prigogine's celebrated principle of minimum, entropy production. Let us recall that, according to 
the latter, physical systems evolve towards stationary states which have minimum entropy production compared to 
slightly displaced neighboring states. Given that the entropy is a convex functional on the state (configuration) space, 
bounded from above for any finite average energy, this implies that the physical evolution will take the system towards 
a local maximum of the entropy or at least towards a ridge. Indeed, in a small enough vicinity of a maximum of 
the entropy, or of a ridge, any evolution with positive entropy production will eventually enter a regime where S 
decreases in time until it vanishes in the equilibrium state or is minimized for the stationary states corresponding to 
a ridge. The variational principle Eq.([T^) only complements this picture by stating that the evolution should follow 
the shortest route to a state of maximum entropy, i.e. the direction of the physical path is selected from among all 
directions satisfying S' > by the requirement that the increase in entropy be maximized at each point in time. In 
this case it can be said that the entropy production evolves towards a minimum of the maximum, to be attained on 
a local maximum or a ridge of the entropy (hyper)surface in state space. 



IV. THE LINEAR NEAR-EQUILIBRIUM LIMIT 



It is natural to anticipate a linear limit for any nonlinear dynamics evolving sufficiently close to a canonical thermal 
equilibrium state, at least in the high-temperature limit. For the modified equation of motion proposed here, the 
linearization process entails essentially the approximation of the entropy operator —plnp to first order in A(p — p"^^) 
around the target equilibrium state 

P^'^^e-P", lnZ = -/3E+^, (42) 
Z ks 

with given average energy E and reciprocal temperature /3. We proceed from the exact expansion 

oo 

-lnp = ^-(/-p)", (43) 

n=l 

which gives for p — p'^'' + Ap, and in symmetrized form, 

- (p^« + Ap)ln(p^« + ^P)^l ^(P^' + ^^)' £ " P'^' " ^'^^"1 ■ 
Separation of the zero- and first-order terms in Ap yields 

- ip'^^ + Ap) Hp'^^ + Ap) = -/^ Inp'^" - 1 {Ap, Inp-^n + ^ ip"', HM} , (45) 

where A(Ap) represents the collection of all terms first-order in Ap from the infinite sum on the right hand side of 
Eq. (p4[) . In order to calculate A(Ap), it is convenient to define the superoperator R and its tilde-conjugate R by 

RAp = (/ - p"«)Ap , (46a) 



RAp = Ap{I - p"") , (46b) 



[R, R] = . 

The expression of A(Ap) can be obtained now in the compact form 
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oo ^ n— 1 

A(Ap) = - 5^ 1 5] R™r"'Ap . 

But R has a well-defined inverse and the super-operator sum in the above expression can be rewritten as 



(47) 



~n-l n-1 



1 m^O 

T OO 



I RR 



oo 



oo /^n 



R R" 

n n 



R R 



E 

n=l 
1 



RR 



In(I-R) -ln(l-R 



Taking also into account that 



and 



where 



p^i • Ap = (I - R)Ap , Ap • = (I - R) Ap 
ln(I - R) = -/3H - (hi Z) I , hi(I - R) = -/3H - (hi Z) I , 
HAp -Ap , HAp = Ap-H, 



we are lead to: 



1 {p^^ A(Ap)} = -| (h - h) . coth [ f (h - H 



Ap 



(48) 
(49) 
(50) 
(51) 

(52) 



Returning to Eq.(^), the first-order in Ap approximation to the entropy operator reads now 

- (p^' + Ap) ln(p^' + Ap) = -p^« In p^« ~ ^{ Ap, H - E} - ^ • Ap - 

z Kb 



(h-h) 



coth 



H H 



Ap 



(53) 



Note that taking the trace in Eq.(53) gives S(t) « S*^' in this regime. Similarly, a simple calculation shows that 
C « (3/2. Assuming also that a « cr'^'' = const. (E, (3) and inserting everything into Eq.(p6[) yields the linearized 
equation of motion 



P = 



H-H 



Ap 



coth 



or , as well, 

A/>=-a^'^^(H-H) 
The general solution of Eq. (^ is given by 

p(t) = e--^'(«*e-°*p(0) + (l - c--^'(«*) p' 



Ap 



where 



^(H-H 



I . 



(54) 



(55) 



(56) 



(57) 



We observe immediately that G is tilde-symmetric, hence it maps any hermitian operator into a hermitian operator, 
and that it preserves probability, since rr[Gp] = 0, rr[e~'^'p(0)] = Tr[p{0)\ = 1. This is entirely sufficient to secure 
the hermiticity of p and the overall conservation of probability. Unfortunately, the action of G does not always 
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preserve positivity and G cannot be identified as a generator of Lindblad type ||T^ . But the positive domain of G 
does include the small neighborhood of p^'' identified as the near-equilibrium domain. Indeed, note first that in the 
diagonal representation of the hamiltonian, the matrix elements of /0°(t) = e^*^'p(0) obey the simple damping law 



pL(t) 



where the (temperature-dependent) relaxation coefficient 7^1, is given by 



|(E^-E,)coth 



1 



(58) 



(59) 



00 

If we consider now an arbitrary state vector |^) — ^ (Eiy|^)|E,y) and the matrix element 



(vl/|p0(t)|vl/) = ^(vI/|E,}p,,(0)(E,|vI/)+ i?e[(vI/|E^}p^,(0)(E,|vI/)] e" 



v=0 



(60) 



it is easily seen that p(t) remains positive for t > if the initial off-diagonal matrix elements p^,y(0), p, ^ v, are 
sufficiently small, as expected for the near-equilibrium regime. On the other hand, one can resort to the equation of 
motion for the state operator 7, Eq(p^), and derive a linear approximation in A7 = 7 — 7'^'^, A7+ — 7+ — (7'^'^)"'' by 
the same procedure as above. The result reads 



A7 = 



H H R 



r) ^ (7"«A7+ + A7(7"«)" 



(61a) 



A7 ■ = - ( a/3 - - ) + 



(h - h) (r - r) ^ (7'='?A7+ + A7 iY")^) 



and shows that, up to first-order terms in A7, 

p(t) = p^i + Ap w [7^9 + A7] (7-^9 )+ + A7^ 



(61b) 



(62) 



such that Ap = A7 (Y'^)^ + Y''Aj+ evolves according to Eq.(|55|) derived above. Furthermore, the conservation of 
energy follows from 



Tr {HAp) = -a'^'^Tr {HAp) 



(63) 



upon recalling that, according to the original equation of motion, the initial state necessarily has the same average 
energy E as the asymptotic equilibrium state. The initial conditions for Eq.(|56|) are so restricted to Tr {HApifi)) = 0, 
which implies of course Tr {Hp{0)) — E. 

As a general feature of the underlying physics, it follows from Eqs.(^6|), ( ^Sf ) and (^9|) that the greater the energy gap 
between two energy eigenstates, the faster the quantum correlation between them is destroyed as the system evolves 
towards equilibrium. On the other hand, the relaxation of the occupation probabilities for each of the energy states 
proceeds at a common rate, independent of the corresponding energy level, since Pvv{^) = ^ j^" ~ ^"'^ ''('^)') -|- 
g-cr ''(/')*p^^(0). As a corrolary, the same holds true for the average of any observable which commutes with the 
hamiltonian, [H, O] = 0, since (0(t)) = Tr[Op{t)] = Tr O e-^^'p(t) e^^' (see Eq.(||)) will involve only p^^-s. 

The same result can be obtained in a formal manner from a generalized Heisenberg representation for Eq.(^), in 
which the observables evolve in time according to 



OA(t) 



a'^iP) G 



/9 



H H 



I 



H H 



OA(t) 



(64a) 



OA(t) = exp 



a^n/3) G 



H-H 



(H-H) 



Oa(0) 



(64b) 
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Here (tl — O ~ [i?, O] and the lower label A reminds that all averages are to be calculated with Ap(0) — p(0)—p'^'^. 
From Eq. (^) above it is immediate that [H, O] = (h. — tlj = yields 

(65a) 



OA(t) -cxp[-a^«t] Oa(0) 



(65b) 



in agreement with the observation above. An unexpected outcome of this result is that the average of an observable 
which commutes with the hamiltonian is conserved throughout the evolution, provided the initial average value is 
identical to the equilibrium average. In other words, {ODeita{0)) implies (0}(t) = {Oy^ . The conservation of energy, 
Eq. (|63[) , is seen to be in fact just a particular realization of this feature. Furthermore, for operators satisfying 
commutation relations of the form [H, A\ = eA, Eqs.(|6^) lead to 



Aa (t) , 



(66a) 



Aa (t) = exp 



(0) , 



where G{x) + 1 = 2; coth(2:). 

Eqs.(|65|) and ( |66| ) allow us to provide a handful of instant examples. 
1) A two- level atom, with the hamiltonian 

i/ = Ei|l)(l|+E2|2)(2| 

and the occupation numbers 

m = (l|p|l} , 712 = (2|p|2} , ni +n2 - 1 , 
obeys a simple relaxation law which follows from Eq.(|65|a): 

ni = -a^«(/3) (ni-n^«(/3)) 

If Eqs.(|67|) are rearranged into the kinctic-like form 

hi = -ki2 ni + k2i n2 , 



(66b) 



(67a) 
(67b) 

(68a) 



h2 = ki2 ni - k2i n2 , 



(68b) 



the corresponding (thermal) transition rates ki2 = f7'^'(/3)n2''(/3) , k2i = a"^"^ {f3)nl'^ {(3) , are seen to have, up to the 
factor of a'^'' , an Arrhenius-like dependence on the temperature. 

2) For a harmonic oscillator of unit mass and frequency lu, described by 

2 2' 

{P)e, = 0, {q)e, = , 
one can apply Eq.(p6|a) to the annihilation and creation operators, 

to recover a coupled system of equations for the average momentum and the average coordinate. 
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(p) = -7(u;, I3){p)~u^q), 



(69a) 



(g) = (p)-7(^, m) 



(69b) 



where 7(w,/3) = cr*^' (/?) [1 + G (/?7iti'/2)]. We recognize a typical damped motion, driven by the classical Langevin 
equation 



{q) + 27 {u, /3) {q) + y + 7' {CJ, /?)] (g) = , 

which is obtained by elimination of the momentum variables from Eqs.(|69|). 
3) For the non-relativistic fr'ee particle hamiltonian 



(70) 



Eg. (pSja) gives the relaxation law 



2m 



(71) 



which shows a (thermal) friction force linear in momentum. In case the initial momentum average coincides with the 
final thermal average, one obviously obtains conservation of the average momentum. More details can be extracted 
from the Wigner function 



/A(p,r,t) = R 
q 

Differentiation of Eq.(|7^) on time and use of Eqs.(p6|), (|5 

/A(p,r-,<) = ^e-^« 

But note that 

^.v,/A(p,f,t) = E«"^'""(4— 



PA(t) 

yields 



P- 



q 



-a^" (/3) ( G ( ) + 1 



i_P_q 
h m 



P-2 



P-2 



PA {t) 



PA (t) 



P.I 



(72) 



(73) 



(74) 



and rewrite the right hand side of Eq.(73) in differential form to obtain 



/a + - • VffA 

m 



_,e,(^)(G(,Mp..V,)+l)/A 



(75) 



The operatorial expression on the right hand side is to be understood in terms of the power expansion G (a;) + 1 = 
X coth (a;) = 1 + 2 ^ (^1)" Cr i'^n) (x / n)'^"' ||l^, where C-r(s) = ^ ^^'^ Riemann zeta function. Hence 

ri=l k=l 



Eq. (|75D reads, in explicit form. 



/a + - • VrfA = (/?) 



n=l 



-l + 2VCi?(2n) (^p- V,- 



2n 



I A 



(76) 



and proves to be a Burnett-type (or generalized Fokker-Plank) equation. Note further that the dependence of /a 
on momentum and coordinate variables can be separated in Eq.([76|), and one can integrate over momentum to 
obtain an exact equation in coordinate space. It can be safely assumed also that the momentum distribution does 
not deviate significantly from equilibrium, such that one can write /a {p,r,t) « nA {f,t) f^"^ {\p\)i where (?% = 
n (r, t) — n^'^ (r, t) is the deviation from the equilibrium value of the localization probability and f^"^ (\p\) is the 
equilibrium momentum distribution. In that case, in the high-temperature limit when only contributions to leading 
order in (3 survive, integration over momentum leads apparently to a diffusion-like equation. 



hA = D(/3) AuA -(T'=«(/3)nA , 



(77) 



12 



with the diffusion coefRcient 



(78) 



where it has been taken into account that Cfl(2) — 7r^/6. But let us recall that Xt — y {h (3)/{i m) is just the de 

Broglie wavelength corresponding to the root-mean-square momentum ^/ip^Y^ , such that in fact D(/3) — a'^'^{P) (At)^- 
It follows necessarily that the diffusion term in Eq.(|7^) can give significant contributions only if the localization 
probability varies substantially on the scale of the thermal de Broglie wavelength At, regardless of the specific value 
of ct'^*(/3). But since states with such variations do not belong to the high-temperature, near-equilibrium regime, we 
are forced to recognize that Eq. ([77|) actually reduces to 



(/?) UA 



(79) 



The linearization procedure developed in this Section can be extended without significant modifications to equilib- 
rium states other than the thermal canonical distribution. It can be shown that the relaxation laws for the elements of 
the density matrix in the diagonal representation of the hamiltonian are similar to those found here for the canonical 
case. A detailed account of this issue will be given elsewhere. 



V. SYMMETRY INVARIANCE, CONSERVATION LAWS AND SEPARABILITY 



It has already been pointed out in Sec. 3 that Eq, 
mations that leave the hamiltonian unchanged. It is a' 



is invariant under any time-independent unitary transfor- 
Iso obviously invariant against time translations, albeit this 
operation can no longer be associated with a unitary transformation. The same is not true, in this form, of time- 
dependent transformations relating different observers in relative motion. But at least in the non-relativistic case, this 
deficiency can be easily corrected so that invariance under the complete dynamical group of the system is recovered. 
Indeed, let us rewrite Eq.(E3) in the form 



-a [plnp + {D (p) ,p}~ pTr (plnp)] + - [p, H] 



(80) 



where D{p) replaces C(p, — E) [H — E), and let us consider the invariance conditions for Eq.(|80|) under a time- 
dependent unitary transformation J7(t), U{t)U^{t) = {t)U {i) = I. As usual, the density matrix becomes 
p' (t) = U{t)p (t) U+ {t), hence 



p'(t) = [/(t)p(t)C/+(t) - [p'(t), ?7(t)C/+(t) , 

while multiplication of Eq. (|80|) by U{t) on the left and U^{t) on the right, followed by use of Eq.(^l]) gives 

p' = -a [p'lnp' + {U{t)D{p)U+{t), p'} - p' Tr (p'lnp')] + 



(81) 



p', U{t)HU+{t)+ihIJ{t)U+{t) 



(82) 



It is easily seen that Eq.(^2|) will regain the form of Eq.(gO|) provided H is invariant under U{t) in the customary 
sense, 



and if, in addition, 



U{t)HU+{t) + ih [/(t)[/+(t) = H 
U{t)D{p)U+{t)^D{p') , 



(83) 
(84a) 



a(p, F - E) = a{p', H-E'), 



(84b) 



In the absence of any evidence to the contrary, the functional a will be assumed in the following to have all necessary 
invariance properties. 
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From Eg. (p3[) it follows in the customary way that if C/(t) spans a Lie group of order n, such that C/(t) = 
exp [{i/h)X^ Kj{t)\ , with A-*, j = 1, 2,. . .n the group parameters and summation over repeated indices understood, 
then the corresponding infinitesimal, hermitian generators Kj{t), j = 1, 2,...n satisfy the familiar commutation 
relations 



(85) 



Note that a conservation law is not yet implied. But let us assume further that the transformations U (t) are such 
that 



U{t)U+{t)^-{a^C,{t) + b) , 



(86a) 



u{t)c,m+{t)^c/ci{t) + f, 



(86b) 



where all parameters a-', b, Cj'' , fj are real functions of the group parameters and time, and the Cj-s are hermitian 
operators (observables). In that case, if the conservation of energy is to be invariant under all transformations U{t), 
it follows from the expression of the transformed average energy 



E 



' = Tr {Hp' (t)) ^E + ihTr(u {t) U+ (t) p' (t)) = E - a^' (Q)' 



(87) 



that a conservation law is required for each Cj. Unfortunately, Eq.(|24D does not account for such supplementary 
constants of motion and simple algebra reveals that D{p) — C,{p, H — Fi) {H — E) does not satisfy the first of Eqs.(p4|), 
despite an invariant hamiltonian, since 

U{t)[C{p,H-E) {H-E)]U+it) = 



C (p, H -ihU (t) u+ (t) ^e) [H -ihii (t) U+ (t) - E 



^C{p\H-E'){H-E') . 

Let us examine now whether modifying Eq. (|8^) to include conservation of the quantities Cj brings about the desired 
invariance under the transformations of the given Lie group. Let the conservation of each Cj be added to the set of 
constraints accounted for in the original variational principle, such that Eq.(n2[) is brought to the form 



6[{j I ln(77+) I 7 ) + ( 7 I ln(77+) | 7 ) + 2 C ( 7 I H | 7 ) + 2 C*( 7 I H | 7 )+ 



(89) 



+ e ( ( 7 I 7 ) + ( 7 I 7 ) ) + 2 r,-'" ( ( 7 |C, I 7 ) + ( 7 I C, |7 ) ) + - ( 7 I 7 ) 







with the new parameters rj^ assumed real, since the corresponding terms will not contribute to the hamiltonian part 
of the equation of motion. Taking again the variation with respect to 7, 7"*" yields 



l7) 



I [In (77+)] I7) + C H|7) + f^^C, I7) + I I7) 



and the corresponding equation of motion for the density matrix, 

p=-<7[p\np+{C{H-E) + {C, - (C,)) , p]~p Tr{p\np)\ + H] . 
Here (Cj) — Tr [Cjp] is the conserved average of Cj, ( and rf are solutions of 

Tr [{H - E)p\np] + 2C Tr \{H - E)^ pi + rj^ Tr [{H - E, Cj - {Cj)} p] = , 



(90) 



(91) 



(92a) 



Tr [{C, - (Q)) p\np] - (t/ha) Tr [ [Q - (Q) ,H~E]p] + 
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+ CTr[{C,-(Q),iJ-E} H+V Tr[{C,-{C,), Ci - {Ci)} p]=0, 



(92b) 



3 = 1, 2,...n 

and we have identified a{ImC) = (l/?i), i?eC ^ C, C = - [Tr (pin p) + 2 (ReC) E + 277-' (Q)] . Eqs.(|9|) always have 
solution, as the matrix of coefficients for the unknowns C, rj^ is recognized to be the positively defined covariance matrix 
for the hamiltonian and the operators Cj . If we presume the invariance of the hamiltonian as defined by Eq. ( ^3| ) , in 
accordance with the discussion above, it is now straightforward to verify that D = ( {H — E) + ij^ {Cj — (Cj)) is 
invariant as well in the sense of Eq.(p4|a), provided (, if change as 

C = C, (93a) 

= 77' c/ + C . (93b) 

In deriving Eqs.(p3|) use has been made of Eqs.(|8^) and the following transformation of {Cj) under the action of U (t): 

(Q) ^ Tr [C,p] = Tr [U {t) C,U+ (t) p'] = c/ (Q)' + /, . (94) 

The complete invariance of Eq.(|9|) requires, of course, that (\ ff defined in Eqs.(p3|) be solutions of the transformed 
Eqs.(p2|), obtained upon substituting p' , E' and {|p^ for p,E and (^P^j respectively. But substitution of p (t) = 
U (t) p' (t) J7+ (t), followed by rearrangement of U, [/+ over observables and use of the relations 

U (t) HU+ {t)^E = H -E' + a^ (Q - (Q)') , (95a) 
U (t) C,U+ it) - (Q) = c/ (Ci - {Ci)) , (95b) 

3 = 1, 2,...7i 
obtained from Eqs.(||), (|§), (|3) and (||), leads to 

Tr [{H - E') p'lnp'] + 2 C Tr [(i/ - E'f p'] + (r/'c/ + ( a') Tr [{H - E' , C, ~ {C,)'] p'] 

+ [i/ha) a^ Tr [[Cj - (Cj)' ,H-E']p']^0, (96a) 
Tr [(Q - (Q)') p'lnp'] - (z/ha) Tr [[C, - (Q)' , H - E'] p'] + 
+C Tr [{C, - {C,y ,H~E'} p'] + {f^"^c,l + (a') Tr [{C, - (Q)' , Q ~ (Q)'} p'] - 

- {t/n<j) a'Tr [[C, - (C,)' , Ci - {Ci)'] p'] = , (96b) 

j = 1, 2,...77 . 

The first of these equations displays the required invariance only if the last term vanishes identically, which demands 

[Q(i),^]-0, (97) 

for 3 = 1, 2, ... 71, while the second equation is seen to be invariant provided 

[Q(t),G(t)] = 0, (98) 

for j = 1, 2,...7i, Z = 1, 2,... 71. We conclude that Eq.(91) is invariant if and only if Eq s.(p| ), ^) and (|9|) 
are simultaneously verified, in which case the parameters C, rj^ transform according to Eqs.(p3^ The generating 
variational principle, Eq. (|89|) , is invariant, of course, under the same conditions. 
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Let us substitute now for f/(t) the special Galilei boost of velocity vq, 



U (t; vo) = exp 



m-X 



(99) 



where m is the total mass of the system, X is the position of the center of mass and P denotes the total momentum. 



Expression (M) obviously prompts the identifications Cj = Pj 



1, b 



(touo)/2, 



mvo 



introduced in Eqs.(|97|) and (|98| ) lead to the recognizable commutation relations 



[H,P,]^[P„Pi]^0. 

Subsequent substitution in Eq.(91) gives the corresponding equation of motion in the form 
p=-a[p\np+{C{H-E)+ rf [P, - {P,)) , p] - p Tr (plnp)] + [p, H] 



which 



(100) 



(101) 



We recover thus in an unexpected manner the celebrated result that the Galilei-invariance of the appropriate non- 
relativistic equation of motion is equivalent to the corresponding invariance of the hamiltonian, the conservation of 
total mom entum and the commutation of the hamiltonian and the total momentum operators. 

Eq.(lOl) reduces to the original Eq.(^J) in the center-of-mass referential, where only states corresponding to an 
eigenstate of zero total momentum for the center-of-mass coordinates need be considered and the dissipative momen- 
tum terms vanish. It also retains the fundamental features previously outlined for Eq.(p4[). In particular, it can be 
checked that pure states evolve according to the usual hamiltonian dynamics, the entropy of mixed states increases 
and the nature of the asymptotic equilibrium states is preserved, up to a slight change of form which accounts for 
the conservation of momentum. It is also evident that Eq.(lOl) is invariant under time-independent symmetry trans- 
formations which leave the hamiltonian and the dissipator D invariant, provided the time-scale parameter a has the 
same property. In particular, if the hamiltonian commutes with the total angular momentum, Eq.( |lOl| ) is invariant 
under finite rotations. However, as for the linear momentum, rotational invariance alone does not imply, in general, 
a conservation law for the angular momentum. The latter can be brought into view by requiring that the equation 
of motion for the density matrix be covariant with respect to all reference frames where the conservation of energy is 
a valid physical law. In particular we should consider translations to observers in uniform rotational motion around 
an axis at rest in some inertial frame. The rather cumbersome details of adding this supplementary constraint will 
be left aside, since nothing new will be gained f or th e formalism. 

A more interesting lack of symmetry for Eq.( |lOl|) , or better, the simpler Eq.(24), lies concealed in the apparent 
absence of separability. Indeed, let the system described by the hamiltonian H be composed of two noninteracting 
subsystems, such that H = Hi + H2, [Hi, H2] — 0, and consider the situation of a separable initial state p{Q) — 
Pi(0)/92(0), of energy E = Ei + E2. Direct inspection of Eq.([2^) shows that the energies of the two subsystems cannot 
be separately conserved and a completely separable solution is thus prohibited. But we also observe that lifting 
the constraint of separate conservation of energy allows a pseudo-separable solution p(t) = pi(t) /32(t) given by the 
coupled system 



Pi 



pi In pi+C <Hi- 



Tr [Hipi] 
Tr (pi) 



Pi 



Pi 



Tr (pi Inpi) 



Tr{pi) 



(102a) 



P2 



P2 In /52 + C i H: 



Tr {H2P2) 
Tr {p2) 



P2 



P2 



Tr{p2 In P2) 



Tr {p2 



^ [P2, H2] 



(102b) 



In this case probability is independently conserved for each subsystem, since Tr(pi) = 0, while energy is only conserved 
globally. 



Tr{Hipi) , Tr{H2P2) 



Tr {pi) 



Tr ip2) 



E . 



The coupling between the (noninteracting) subsystems appears to be as instantaneous and nonlocal as usual quantum 
entanglement, but unlike the latter, it involves an unorthodox exchange of energy. The significance of this unusual 
outcome follows from the observation that, according to Eqs.(102), the equilibrium of the compound system is attained 
for values of a and ^ common to both subsystems, hence for a common generalized temperature. Imagine now that 
the initial states for the two subsystems are chosen as individual equilibrium states with different corresponding 
temperatures. It follows that the dynamics given by Eq.(|23) will drive the total system towards a new state of 
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equilibrium, with a temperature common to both components. We cannot but concede the obvious similarity of 
this unconventional effect with the classical process of equilibration by thermal contact. Its origin lies in the very 
assumption of maximal entropy increase on which Eq.(p^) has been derived. Indeed, even when the entropy of each 
subsystem is already maximal under individual isolation, if states of larger total entropy are available, probabilities 
and energy (heat) will be necessarily redistributed so as to enforce a further increase of the overall entropy. Whether 
this entropic entanglement, or ideal thermal contact, is or not an element of reality appears equivalent to accepting 
or rejecting the conjecture that an isolated, perfectly ideal gas can undergo relaxation towards equilibrium. 

We can provide formal support towards the positive by pointing out that the effect of entropic entanglement does 
not necessarily interfere with the concept of separable evolution for mutually isolated systems. First let us note that 
explicitly specifying an adiabatic separation (in the thermodynamic sense) of the noninteracting systems, and hence 
allowing for separate conservation of energy, removes most of the entropic entanglement. In this case the resulting 
equation of motion will display distinct ^-s for each of the systems, but a common time-scale parameter, i.e. 



with 



p = -CT [p In p + Ci {Hi -Ei,p} + C2 {H2 -E2,p]~pTr{p In p)] + - [p, + H2] , (103) 



Tr(H,p) 
Ei = m ~ const. , 



Trip) 



for i — 1,2. As before, it proves possible to extract a pseudo-separable solution p{t) = pi{t)p2(t), but Eqs.(102) are 
replaced by 



Pi = -c 



pi Inpi + Ci {Hi - El , pi} - pi 



Tr (pi Inpi) 
Tr{pi) 



[pu Hi] , 



(104a) 



+ ^ [P2, H2] . (104b) 

where this time the parameters, i — 1,2, will be found to depend only on the corresponding pi and Hi, in exactly 
the manner obtained for a single isolated system. Yet the two evolutions remain tethered by the time-scale parameter 
cr, thus retaining a weaker form of entropic entanglement. The simple presence of other noninteracting, adiabatically 
separated systems appears to alter the time-scale of dissipative relaxation for any given system. If a is assumed 
variable in time, e.g. through a dependence on p, this influence will be time-dependent unless all other systems have 
reached equilibrium. But since a does not affect the nature of the asymptotic equilibrium state, the equilibrium of 
any one system will not be disturbed by other systems and will display an individual temperature determined solely 
by the corresponding energy content. 

A careful examination will trace the above type of nonseparability to the fact that the corresponding variational 
principle selects the direction of maximum entropy increase by referring to the time derivative of the total (entangled) 
state operator, and not to disentangled, individual state operators separately. However, this pitfall can be avoided if 
it is recognized that true mutual isolation precludes entanglement on invariance grounds. Indeed, regardless of the 
nature of the underlying dynamics, the evolution of two mutually isolated systems should remain invariant under every 
transformation pertaining to the individual symmetry groups. In particular, it should be invariant under individual 
time translations. Since entangled states certainly do not possess this invariance, they do not describe truly isolated 
systems. In other words, the restricted subspace of the state space that can be spanned by the dynamics of mutually 
isolated systems should contain only non-entangled states and the evolution of each of the factor states should be 
driven independently. In our nonlinear setting, where this subspace is selected by means of the generating variational 
principle, this restriction has to be correctly built in the variational functional itself. Hence one has to account both for 
individual conservation laws, excluding thus any energy exchange, as well as for vanishing entanglement. The latter 
imposes a separable state operator 7(f) = 7i(t)72(i) and also requires that the entropy production be maximized 
separately with respect to variations of 71 and 72, i.e. the a term in the variational principle should be replaced 
according to 

2 2 2 

-(7I7) — (71 72 1 71 72 ) H (7172I7172) 

cr cri cr2 

with each cr^ a functional only of 7^ and Hi. But then the variational principle takes the form 

'5{(72|72)J^i + (7i|7i)^2} = , (105a) 



P2 = — cr 



P2 In p2 -I- C2 {H2 - E2 , P2} - P2 



Tr (p2 In P2) 
Trip2) _ 
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= (7»| ln(7^7+)l7.) + h^\H7^lt)\i^) + 2 C^{j^mh^) + 2 C;(7»|H.l7.)- 



(105b) 



where 



i = 1, 2 , 



6 = 6 + U2 + C2E2 - 



ke (72 1 72 J 



Si 



kB(7i|7i^ 



Independent variation on 71 and 72, followed by extraction of the Lagrange parameters from the corresponding 
conservation conditions leads now to the desired separate equations of motion for pi — "fi"ff , 



Pi In Pi + C,i{Hi - E,, Pi} - p., 



Tr{pi In Pi) 



Tr{p^) 



^ [Pi. Hz] 



(106) 



I = 1, 2 , 

with 

CTi = (7i{pi,Hi) > , 

_1 rr[(g, -Ei)pilnpi] 
2 Tr[{H,-Ei)^pi] ' 

rr(g,p,) 
Ei = — — — ^ = const. . 

Obviously, the invariance of the nonlinear dynamics under the symmetry group of each component subsystem is so 
restored, provided the tii-s are invariant also. 

VI. GENERALIZATION TO ARBITRARY ENTROPY AND ENERGY FUNCTIONAL FORMS 

The framework developed in the previous Sees, can be easily expanded to accommodate non-standard entropy 
functionals and/or energy forms with a nonlinear dependence on the density matrix p. This generalized formalism can 
then provide nonlinear extensions for, e.g., the Lie-Poisson dynamics or a standard hamiltonian evolution supplemented 
by a nonextensive Tsallis entropy appropriate for systems with fractal properties. We sketch here only the 
derivation of the generalized equation of motion, since a detailed analysis exceeds the purpose of the present work. 

To this end, let us start with a Lic-Poisson equation of motion in the form |0] 



P = 



P, H{p) 



(107) 



where H{p) is in general a hermitian, nonlinear functional of p. The energy conservation law is now replaced by 

Tr (h{p)p) = 

or in terms of the state operator 7, 



(7|H(p)l7) + (7|H(p)|7) = 



(108) 
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The law of probability conservation, on the other hand, remains unchanged since Tr{p) = or 

(7I7) + (7I7) = 



(109) 



Let us search now for a nonhnear evolution that observes the above conservation constraints, Eqs.( |lOq ) and (109) 
and is also subject to a second principle based on some unspecified, positive definite entropy functional = TrlSlp)) 

Kb 

such that S = Tr {{5SI5p)p \ > or 



/ . , (^S , , , , JS , . , 

(7l^l7) + (7l^l7) >0 
op op 



(110) 



Here the operator S{p) is assumed hermitian and {SS/6p) denotes its hermitian functional derivative with respect to 
p. The corresponding variational principle is now written 

S |-(7 I^l7) - (7l^l7) + 2 C (7|H(P)|7) + 2 C (7|H(p)l7) + C [ (7I7) + (7I7) ] + ^ (7I7) | = . (Ill) 
and can be verified to generate the following equation of motion: 



P - 



-^p + C {Hip)-{H{p)), p} 



'dS\ 



h 



P, H{p) 



(112) 



where 



and 



(A) 



Tr{Ap) 



Trip) 

[Hip)-{Hip)))iSS/6p) 



[Hip) ^i/(p); 
a = aip,Hip)-(Hip))) >0 



We note that if i5S/Sp)p = for pure states, p = p^ , then (SS/Sp) = 0, C = 0, and the pure state dynamics reduces 
to that prescribed by Eq.( |l07] ). 

When the energy functional reduces to the hamiltonian, Hip) — H, and the entropy is given the standard von 
Neumann expression, such that Sip) = — plnp, idS/Sp)p = —p\np — p, we recover the basic Eq. (p^) . A p-dependent 
Hip), complemented by the standard entropy, leads to a nonlinear extension of the Lic-Poisson dynamics, 



p\np + C [Hip) - (^Hip)) , p} 



Tr (pin p) 
Trip) ' 



h 



P, Hip) 



(113) 



with 



Tr 




[H-{Hip))) 


plnp 


Tr 




[Hip) (Hip))) 


2 

P 



C = -: 



If Hip) is reduced to the standard hamiltonian H, but the entropy is given a Tsallis form, with 



Sip) 



P-P" 

g-1 
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5S p- qp'i 
dp q~l 

for given real q, the result will be a nonlinear extension of the von Neumann dynamics under Tsallis q-thermostatistics, 
which reads, after a few elementary manipulations, 



' P^ + C{H-E,p}-^^p 



where 



q-l ' g - 1 Tr{p} 

1 q Tr [{H - E)^^] 



J [P, H] , (114) 



C 



2q-lTr[{H - E)2p] 



Situations where the standard averages have to be replaced by q-averages can be approached in the same fashion, by 
appropriately redefining the conserved functionals. 

VII. CONCLUSION 

We have constructed and analyzed a non-relativistic nonlinear extension of the quantum law of evolution, which 
accounts for the second principle of thermodynamics and is not at odds with the factual linearity of pure state 
propagation. The theoretical existence of such an extension confirms that the linear and unitary evolution of pure 
states is not in itself sufficient proof for the general linearity of quantum mechanics H,^. One must conclude that 
the linear propagation of mixed states also has to be corroborated experimentally, to comparable precision, before a 
definitive conclusion can be drawn. It is hoped that the formal study developed here provides a meaningful benchmark 
in this sense. 

Our main result is Eq. (p^ , which defines the modified time evolution of the density matrix. The equation of 
motion was extracted from a variational principle on the space of state operators, rather than the space of density 
matrices, as a trajectory of maximal entropy production under the constraint of energy and probability conservation. 



augmented eventually by the requirement of Galilei invariance (see Eq.(lOl)). Should we drop the requirement of 
entropy increase, the parameters -ReC, C vanish and the equation of motion reduces automatically to the common 
hamiltonian form. The outlined procedure may not be unique, but is encouraging in its consistency. In addition, it 
applies as well to alternate theories which use nonstandard energy or entropy forms. It is notable that the variational 
principle has sense only in terms of state operators, whereas the equation of motion can be stated simply in terms of 
the conventional density matrix. 

A peculiar and unexpected idea brought forth in our ansatz is that a maximal increase of entropy does not necessarily 
result in maximal decoherence, to the effect that a pure state of a perfectly isolated system is not allowed to evolve into 
a mixed state. On the contrary, the proposed quantum equivalent of the second principle of thermodynamics is seen 
to introduce only a limited degree of decoherence, in the sense that the cardinality of the set of nonzero eigenvalues 
of the density matrix is preserved. As already mentioned, for the particular case of a pure initial state this leads to 
the usual unitary evolution. The same property also supports, aside from canonical equilibrium states, a rich class of 
"negative-temperature" equilibrium states, which bring to mind the notion of thermal coherence. Furthermore, the 



ideal thermal contact phenomenon discussed in Sec. 5 abides by the same rule and, according to Eqs.(102), a system 
in an initially pure state will remain in a pure state even if it is in contact with, but not necessarily interacting with, 
other systems. However, in that case the pure state und ergoes relaxation according to a dynamics of Gisin type [Q, 
as seen by taking, e.g., pi — p\, pi Inpi = 0, in Eq.(|l02|a), 



P. = -.c{i/i-^^;^,Pr} + J[Pr,i/.]. (115) 

Depending on the sign of the asymptotic stationary state is an energy eigenstate for the lowest (if ^ > 0) or for 
the highest (if C < 0) energy level contributing to the initial state pi(0). As detailed in Sec. 5, the state of thermal 
contact is not to be mistaken for a state of mutual isolation, despite the absence of explicit interactions. 

We find it promising that bending quantum dynamics to account for classical phenomenological irreversibility 
suggests a rather unified picture of both reversibility and irreversibility, as well as coherence and decoherence, while 
preserving such fundamental features as symmetry invariance. However, the self-consistency of the theory is limited 
at this point by the need for an explicit expression for the entropy production, which means that the equation of 
motion remains determined up to the scale setting functional a. We leave the resolution of this problem for future 
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consideration, although a definite expression for <j certainly conditions the consistency of our results. For instance, 
Eq. ( ^9|) for the near-equilibrium damping constants of the density matrix elements between energy eigenstates shows 
an acceptable dependence on the energy gap between the states, but the wrong temperature dependence (7^,^ — > 
as /3 — > , 00 as P —>■ 00) if a is assumed temperature independent. In the least, this observation serves to 

hint that a should behave like f3~^'^~^^\ (5 > 0, in the vicinity of canonical equilibrium, which in turn can be used, of 
course, as a theoretical benchmark. 
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